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ABSTRACT 

We present a current best estimate of the integrated near-infrared (NIR) extragalactic background 
light (EBL) attributable to resolved galaxies in J, H, and K s . Our results for measurements of vlu 

in units of nW m~ 2 sr _1 are 11.7^2 g in J ■, H-Sli^ m H anc ^ m -^s- We derive these new 

limits by combining our deep wide-field NIR photometry from five widely separated fields with other 
studies from the literature to create a galaxy counts sample that is highly complete and has good 
counting statistics out to JHK S ~ 27 — 28. As part of this effort we present new ultradeep K s — band 
galaxy counts from 22 hours of observations with the Multi Object Infrared Camera and Spectrograph 
(MOIRCS) instrument on the Subaru Telescope. We use this MOIRCS K s — band mosaic to estimate 
the total missing flux from sources beyond our detection limits. Our new limits to the NIR EBL 
are in basic agreement with, but 10 — 20% higher than previous estimates, bringing them into better 
agreement with estimates of the total NIR EBL (resolved + unresolved sources) obtained from TeV 
7— ray opacity measurements and recent direct measurements of the total NIR EBL. We examine 
field to field variations in our photometry to show that the integrated light from galaxies is isotropic 
to within uncertainties, consistent with the expected large-scale isotropy of the EBL. Our data also 
allow for a robust estimate of the NIR light from Galactic stars, which we find to be 14.7 ± 2.4 in J, 
10.1 ± 1.9 in H and 7.6 ± 1.8 in K, in units of nW m" 2 sr -1 . 



Subject headings: cosmology: observations 
background) 



galaxies: fundamental parameters (counts, near-infrared 



1. INTRODUCTION 

The near-infrared (NIR) extragalactic background 
light (EBL) is the total light from resolved and unre- 
solved extragalactic sources in the NIR. This represents 
the integrated light from all star and galaxy formation 
processes over the history of the universe that has been 
rcdshifted into the NIR. Some fraction of the NIR EBL 
can be resolved as the light from individual galaxies (In- 
tegrated Galaxy Light, IGL), but the fractional contribu- 
tion from unresolved (and perhaps unresolvable) sources 
is not well constrained. A measurement of the total NIR 
EBL minus the IGL will provide insight into the energy 
budget of the early universe. 

While direct unresolved measurement of the NIR 
EBL is technically difficult due to complex fore- 
grounds, several authors have reported a measure- 
ment of the total (resolved plus unresolved) NIR EBL 
(Mauser et al. Ill 9981: [Dwek fc Arendtlll998HGoriian et alJ 
' 20001: iWrightlboOll: iMatsumotol 120011: IMatsumoto et al.l 
2003 ICambresv et al.l 120011: iLevenson et al.l I2007D and 



found it to be a factor of two o r more above the IGL ob- 
tained throu gh so urce counts dMadau fc PozzettH 120001 : 
iTotani et afl 120011: iThompsonl 120031 ) . This is known as 
the NIR background excess (NIRBE). The spectral en- 
ergy distribution (SED) of this measured excess, in some 
cases, appears very similar to that of zodiacal light, which 
suggests there may be a foreground contamination is- 
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sue. Another possible solution to this excess is a large 
population of undetected faint galaxies and/or popula- 
tion III (PopIII) stars contributing a large fraction of the 
NIR background (see [Kashlinsky 2005; Haus er fc Dwekl 
2001 for reviews). For example. IMatsumoto et al.l (|2005l ) 
found a large NIRBE with a distinct spectral break at 
~ 1/im, which they interpreted as a possible signature 
of a Lyman limit break in PopIII stellar spectra at a 
redshift of ~ 10. 

Since the NIR EBL presents a source of opacity for 
TeV 7— rays via pair production, the density of the back- 
ground light can, in principle, be measured via direct ob- 
servation of TeV blazars. This method relies on assump- 
tions of the intrinsic blazar spectrum and the SED of 
the EBL from ultraviolet to the NIR, both of which are 
poorly observationally constrained. Nevertheless, TeV 
7— rays provide an independent estimate of the NIR EBL 
that can be used to help determine what fraction of the 
background light can be attributed to resolved sources, 
and how much may arise from faint and possibly exotic 
s ources in the early u niverse. 

iDwek et alJ (|2005|) explored the possibility that the 
NIRBE is of extragalactic origin by looking for the ab- 
sorption imprint of the ~ 1/im spectral break on the 
7— ray SEDs of several TeV blazars. They concluded 
that the apparent NIR excess is likely not of extra- 
galactic origin because its signature is not detected in 
the blazar SEDs. Recent observations of TeV blazars 
with the High Energy Stereoscopic System (HESS) have 
placed limits on the NIR EBL that allow for the possibil- 
ity of little or no exces s above integrated galaxy counts 
(lAharonian et al 1 12006). 

iMazin fc Raud (|2007l ) used 13 TeV blazars and a grid 
of NIR background SEDs to further constrain the back- 
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ground light. They derived limits on t he NIR EBL that 
are ro ughly in agreement with those of lAharonian et all 
(2006). However, systematic uncertainties in these 
7— ray opacity measurements due to uncertainties in the 
blazar SED and that of the EBL result in large error 
bars that could still allow for a significant NIRBE (of 
the same order as the IGL) to be present, even if their 
assumptions about the intrinsic blazar SED and that of 
the EBL are assumed to be reasonable. 

As such, the goal of this paper is to derive a firm 
lower limit to the NIR IGL so as to facilitate a ro- 
bust comparison with future new results from direct 
and indirect measurements of the NIR EBL. We present 
the best current estimate of the NIR IGL by combin- 
ing our deep wide-field NIR data set with a compila- 
tion of galaxy counts from the literature and ultradeep 
K s — band galaxy counts from a subfield of the Su baru 
MOIRCS image p resented in iBarger et all 1)20081 ) and 
IWang et al.l (120091) . 

The structure of this paper is as follows: We describe 
our NIR dataset briefly in Section [2j We calculate the 
contribution of resolved galaxies to the NIR EBL in Sec- 
tion [3] We compare with previous measurements of the 
NIR EBL and IGL and discuss the implications of our 
results in Sectional We summarize in Section [3] Unless 
otherwise noted, all magnitudes given in this paper are 
in the AB magnitude system (toab = — 2.5 log /„ —48.60 
with /„ in units of ergs cm~ 2 s _1 Hz -1 ). 

2. THE NIR DATA 

We presented the observations, data reduction, star- 
galaxy separation, and b right galaxy counts f or our deep 
wide-field NIR dataset in lKeenan et all (|2010H (K10 here- 
after). The data consist of five widely separated fields 
at high galactic latitudes, including one centered on an 
Abell cluster (Abell 370). The survey covers approx- 
imately 3 deg 2 reaching a 5 a limiting magnitude of 
JHK S - 22 - 23 over - 2.75 deg 2 with another - 0.25 
deg 2 to JHK S ~ 24, making it one of the deepest wide- 
field surveys to date in the NIR. The fields covered in 
this survey are described in Section |2"7T1 

In this paper we supplement the above data with those 
from an ultradeep K s — band mosaic composed of imaging 
done by our group and Japanese groups led by various 
investigators using the Multi Object InfraRed Camera 
and Spectrograph (MOIRCS) instrument on the 8.2 m 
Subaru Telescope. The MOIRCS imaging is described in 
Section O 

2.1. Wide-Field Imaging and Galaxy Counts 

Our first two fields are centered on t he Chandra Large 
Area Synoptic X-ra y Survey (CLASXS: lYang et al.lf2003: 
iSteffen et al.l [2004s) an d the Chandra Lockman Area 
North Survey (CLANS: iTrouille et al.ll2008L [2009?) . Each 
of these fields cover ~ 1 deg 2 in JHK S . These fields 
are located in the Lockman Hol e region of extremel y 
low Galactic HI column density ([Lockman et al.|[l986l) . 
Our third field covers a 0.25 deg 2 a rea centered on the 
Chandra Deep Field N orth (CDF-N, iBrandt et aLll2001t 
I Alexander et al.ll2003D . The CDF-N contains the Great 
Observatories Origins Deep Survey North (GOODS-N; 
145 arc min 2 HST Advanced C amera for Surveys obser- 
vation, Giava hsco et al.ll2004D . Our fourth field is the 



Abell 370 (A370) cluster and surrounding area (~ 0.5 
deg 2 ). A3 70 is a cluster of richness at a redshift of 
z = 0.37. Our fifth field is a ~ 0.2 deg 2 area centered 
on the "Small-Survey- Ar ea 13" (SSA13) fro m the Hawaii 
Deep Fields described in iLillv et all ([19911 ). 

In K10 we used the bright NIR galaxy counts 
from our survey in combination with counts we 
generated from The Two M icron All Sky Survey 
(2MASS, ISkrutskie et "all [20061). The UKIRT Inf rared 
Deep Sky Survey (UKIDSST lLawrence et all 12001 and 
other selections from the literature to explore local large 
scale structure via the slope of the galaxy counts curve 
as a function of position on the sky. Here we extend the 
counts from K10 to fainter magnitudes using the meth- 
ods described below. 

We generated average galaxy counts for our ensemble 
of NIR images by the methods described in K10. We cor- 
rected the counts for completeness through a simulation 
in which relatively bright (JHK S = 19) galaxies were ex- 
tracted from each image and then scattered at random 
to other parts of the image to test how well they were re- 
covered as a function of apparent magnitude. We faded 
the test galaxies in 0.25 magnitude steps until the recov- 
ered fraction fell below 25%. We only included counts 
in the final average where completeness was better than 
50%. 

In our completeness simulations, given that we know 
the total magnitude of objects we are scattering through- 
out the image, we are able to compare the recovered 
magnitudes in various apertures to total magnitudes and 
thereby derive a magnitude correction factor for faint 
objects. Figure [T] shows example results from this ex- 
ercise. In Figure QJ, we show the fractional complete- 
ness (N recovered / 'N scattered) as a function of apparent 
magnitude for the CDF-N K s — band. In all fields we 
reliably recover, on average, relatively bright objects 
(JHK S < 20) at the correct magnitude using the SEx- 
tractor MAG_AUTO, which fits a Kron-like elliptical 
aperture to each object ( Kron 1980). We find that the 3" 
aperture shows a roughly constant offset (~ —0.2) from 
total magnitude. Based on these results, we measure to- 
tal magnitudes in the MAG_AUTO aperture for objects 
brighter than 20 th magnitude and switch to corrected 3 
arcsecond aperture magnitudes for objects fainter than 
20 t/l magnitude using the measured difference at 20 th 
magnitude. This magnitude correction only affects our 
calculated IGL at the 0.1 nW m -2 sr _1 level, which is 
roughly a factor of 10 below the uncertainty in the IGL 
due to the uncertainty in the counts themselves. 

2.2. MOIRCS Imaging and Galaxy Counts 

The Subaru MOI RCS K R — band image i s centered on 
the GOODS-N field ([Giavahsco et al.ll200l. Part of the 
Japanese data were published in lKaiisawa eFall ((2006h . 
Here we include all data taken between 2005 and 2008. 
T he data reduc t ion fo r the MOIRCS mosaic is described 
in IBarger etHI ([2008D and IWang eFall jjp)!): and is es- 
sentially the same as that used in K10. 

K s — band galaxy counts from the MOIRCS mosaic 
presented here are taken from a subfield of ~ 15 
arcmin 2 with an integration time of ~ 22 hr yielding 
a 5 a limit of ~ 26.5. This MOIRCS image is, there- 
fore, approximately as deep as the Subaru Deep Field 
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Fig. 1. — Example results from our completeness simulations 
(this case is for the CDF-N K s — band, though results are similar 
across all fields), (a) Fractional completeness as a function of ap- 
parent magnitude. Asterisks and solid black line show fractional 
completeness (N recoveret i/N st . atttirtil i) as a function of apparent 
magnitude (b) Magnitude difference as a function of apparent mag- 
nitude comparing the MAG-AUTO and 3" apertures and showing 
total magnitude minus recovered magnitude for the two different 
apertures. Black squares show results for the MAG_AUTO aper- 
ture and blue diamonds for a 3" aperture. The MAG-AUTO aper- 
ture does a good job at recovering total magnitude for objects 
brighter than 20 th magnitude, and the 3" aperture shows a roughly 
constant offset (~ —0.2) from total magnitude. 

(SDF. lMaihara et al.l200lD a nd Subaru Super Deep Field 
fSSDF. IMinowa et alJ 120091 but more than an order of 
magnitude larger in area on the sky than either survey. 

We generated a catalog for the MOIRCS field us- 
ing SExtractor software version 2.4.4 (jBertin fc Arnoutsl 
Il996f ) for source identification and photometry. We used 
the MAG_AUTO aperture and set the minimum number 
of contiguous pixels to be considered a detection (SEx- 
tractor parameter DETECT-MINAREA) to 2 in order 
to push the counts on the faint end. We did not remove 
stars from the MOIRCS counts because we begin count- 
ing galaxies at 23 rd magnitude where stars identifiable 
through color selection do not contribute significantly to 
galaxy counts (~ 1% for JHK S = 23 from K10 star 
counts compared with average galaxy counts from this 
paper) . Cool dwarf stars unidentifiable through color se- 
lection may contaminate the counts on a slightly higher 
level, but their space density is unknown and we have no 
means of separating them from galaxies. 

For the Subaru MOIRCS K s — band data, we do not 
correct for completeness or total magnitudes. Rather, 
we calculate the total missing flux in the image using 
a met hod similar to that employed by iThompson et al.l 
((2001 . We discuss this in more detail in Section 13.11 
We adopted this method because it allows us to not only 
account for the flux in objects missed near our detec- 
tion limit (completeness correction) , but also for the flux 
missed in objects or parts of objects below our detection 
limit. In J— ba nd and H— band we ad opt the missing 
flux estimates of IThompson et all (|2007D . 



3. NIR BACKGROUND DUE TO RESOLVED GALAXIES 

In Figure [5] we show our completeness-corrected and 
averaged galaxy counts with those drawn from the lit- 
erature. However, before combining all of our data in 
this way, we first investigated field to field variations 
for the IGL for our five fields over the magnitude range 
14.5 < JHK S < 22.5 where all five are highly complete 
and have good counting statistics. We found the IGL 
over this magnitude range to be consistent across the 
four non-cluster fields (CLASXS, CLANS, CDF-N and 
SSA13 - 7 - 8 nW m~ 2 sr -1 ) with altj dispersion of 
±0.5 nW m~ 2 sr -1 . As such, we find that the IGL is con- 
sistent with l arge-scale istropy , an expected signature of 
the EBL (see lKashlinsky||2005l: lHauser fc Dwekll2001l for 
reviews). 

In the A370 cluster field, our IGL measurements were 
^2 — 3 nW m~ 2 sr -1 higher in all bands with the peak 
contribution to the excess light arising from galaxies at 
JHK S ~ 17, consistent with an excess of L* galaxies at a 
rcdshift of ~ 0.4 as in A370. We include the cluster field 
in our average counts for this study because in a survey 
of a few square degrees such as this, roughly one massive 
cluster will be present. 

In Figure [2] we use black diamonds to denote the counts 
for our five wide fields and black squares to denote counts 
for the MOIRCS image. The counts determined by K10 
from the 2MASS field with Galactic latitudes of |6| > 30 
are denoted by red asterisks. The counts determined 
by K10 from the the 2MASS-6x Lockman Hole survey 
are denoted by green triangles. The counts determined 
by K10 from three subfields of the UKIDSS LAS are 
denoted by blue plus symbols. The three subfields that 
make up this section of the LAS are described in detail in 
K10. Other data points are taken from the publications 
listed in the plots. The dashed curves represent the error- 
weighted least squares running average from which we 
calcuate the IGL. We describe the methods we used to 
calculate this average in Section 3 

Figures EJd-f) display the same data as in panels (a- 
c) after dividing by an arbitrarily normalized Euclidean 
model (of a = 0.4 in the form N(m) — Ax 10 am , and A 
is a constant) to expand the ordinate and demonstrate 
where resolved galaxies contribute the most to the IGL. 
A flat line in (d-f) would imply an equal contribution to 
the IGL at all magnitudes. The areas of positive slope 
show where galaxies contribute a larger fraction to the 
IGL as one moves toward fainter magnitude. The steep 
negative slope beyond JHK S > 23 — 24 demonstrates the 
diminishing contribution of resolved galaxies to the IGL 
at the faintest magnitudes. 

To calculate the IGL from all the counts data displayed 
in Figure[2j we used an error-weighted least-squares fit of 
the the slope and normalization parameters "A" and "a" 
in the function N(m) = A x 10 am . We fit the data over 
a range of ±2 magnitudes above and below each magni- 
tude bin to get two estimates for the counts in that bin. 
In other words, at a given apparent magnitude "m" , we 
fit the range m — 2 to m and m to m + 2 separately. 
Each fit yielded an estimate of the counts at apparent 
magnitude m. We then recorded the average of these 
two values as the best estimate for counts in that bin. 
Each of the two fits also yielded a 1 a confidence interval 
for the counts in each bin. Given these fitted values and 
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Fig. 2. — (a-c) Galaxy counts as a function of apparent magnitude. Average galaxy counts from our deep, wide-field data on 5 fields (this 
study) are denoted by black diamonds. The K s — band counts from our Subaru MOIRCS data (this study MOIRCS) are denoted by black 
squares. Error bars for this work are approximately the size of the plot symbols. The counts determined by K10 from the 2MASS field with 
Galactic latitudes of |fe| > 30 are denoted by red asterisks. The counts determined by K10 from the the 2MASS-6x Lockman Hole survey 
arc denoted by green triangles. The counts determined by K10 from three subfields of the UKIDSS LAS are denoted by blue plus symbols. 
Other data points are taken from the studies listed in the plots. The dashed curve shows our error- weighted least squares running average 
(described in Section Q5{ from which we calculate the NIR IGL. (d-f) The same data as in (a-c) but divided by an arbitrarily normalized 
Euclidean model with slope a = 0.4. 



confidence intervals we found the maximum and mini- 
mum allowable counts over the 1 a range and recorded 

maximum-minimum &s ^ fin&1 j ff errQr associatec l wit h 

the average value for the counts in that bin. We then 
integrated these fitted values for the galaxy counts to 
calculate the IGL and its associated 1 a confidence in- 
terval for the magnitude range 11 < JHK S < 28. Our 
fitted values for the galaxy counts with the confidence 
intervals as described above are given in tabular form in 
Table [U 

The results of this calculation are shown in Figure |3l In 
Figure EK, the solid curves represent the cumulative con- 
tribution of galaxies to the IGL as a function of apparent 
magnitude, and the dashed curves show the percent con- 
tribution in each one-magnitude bin. Figure [3Jd shows 
the calculated slope of the galaxy counts curve as a func- 
tion of apparent magnitude. Due to simple geometrical 
effects, when the slope drops below a = 0.4 the total 
light from galaxies begins to converge. 

Figure [3b shows that the counts slope for all three 
bands begins near the expected Euclidean value of a = 



0.6 at bright magnitudes and then drops steadily beyond 
JHK S ~ 17. Due to these trends in slope, the vast ma- 
jority (~ 90%) of the resolved IGL arises from galaxies in 
the range ~ 15 < JHK S < 24. Over this entire range our 
study is highly complete and has good counting statis- 
tics. 

Furthermore, from Figure [3] it can be seen that for 
sources fainter than JHK S = 28 to make any apprecia- 
ble contribution to the IGL the counts curve would have 
to rise dramatically over several magnitudes beyond the 
limits of current surveys. Such a rise would quickly re- 
sult in several (or more) galaxies per square arcsecond, 
rendering any counting exercise impossible due to confu- 
sion. Faint galaxies certainly exist beyond the limits of 
the deepest NIR surveys, because at high redshifts the 
faintest apparent magnitudes observed are only probing 
a few magnitudes fainter than L» down the luminos- 
ity function. It is unknown whether there is a steep 
upturn in the luminosity function toward faint magni- 
tudes, but if such an upturn exists and faint galaxies 
contribute significantly to the NIR EBL, they would need 
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TABLE 1 

Average Galaxy Counts Used to Calculate the NIR EBL 
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1.851 
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3.340 
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2.644 
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3.117 


3.825 


2.736 


3.901 


2.479 


21 


4.098 


3.495 


4.148 


3.368 


4.200 


2.701 


22 


4.455 


3.811 


4.488 


3.476 


4.459 


2.946 
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4.719 


4.162 


4.781 


4.059 


4.795 


3.669 
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4.916 


4.391 


4.965 


4.170 


5.008 


3.931 
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5.067 
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4.732 


5.178 


4.653 


26 


5.285 


4.439 


5.297 


4.647 


5.324 


4.618 


27 


5.370 


4.297 


5.361 


5.055 


5.462 


5.227 


28 


5.376 


4.717 


5.456 


4.755 


5.547 


4.983 



a AB Magnitude. 

b N is the surface density of galaxies in units of mag _1 deg~ 2 . 

c AN is the 1 a error estimate for N (in the same units) derived from the fits described in 
Section [3] to all the data shown in Figure [2] 




10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 
AB Mognitude 



Fig. 3. — (a) Percent contribution to the IGL as a function of ap- 
parent magnitude. Solid curves show the cumulative contribution 
from galaxies to the NIR IGL as a function of apparent magni- 
tude. Dashed curves show the fractional contribution per each 
magnitude bin. This demonstrates that the vast majority of the 
IGL (~ 90%) arises from galaxies in the apparent magnitude range 
~ 15 < JHK S < 24 where this study is highly complete and has 
good counting statistics, (b) Slope of the galaxy counts curve as 
a function of apparent magnitude. Dash-dotted curves show the 
measured slope of the galaxy counts curve as a function of appar- 
ent magnitude. When the slope drops below a = 0.4 the total light 
from galaxies becomes convergent. 



to be so numerous as to be unresolvable. As such, it 
may be safe to say that the resolvable portion of the 
NIR EBL (the IGL) has, for the most part, been mea- 
sured and that the most important contribution to the 
resolved portion comes from galaxies in the magnitude 
range ~ 15 < JHK S < 24, for which the deep wide-field 
data presented here are optimized. 

3.1. K s — band Faint Counts Slope and Missing Flux 



We use our ultradeep MOIRCS data to test for a pos- 
sible steepening of the counts slope at faint magnitudes 
in the K s — band, as appeared po ssible from the SDF 
(|Maihara et all l200l and SSDF ()Minowa et al.l l2fJ05l) 
K s — band counts. We see from Figure^, that the slope 
of the counts in the MOIRCS field continues to become 
natter out to the faintest magnitudes. We detect no up- 
turn in the slope for very faint galaxies. As mentioned 
above, the MOIRCS counts have not had a complete- 
ness correction applied. Such a correction would steepen 
the curve slightly at the faintest magnitudes, but this ef- 
fect would not be strong enough to alter significantly the 
tr end in slope. 

iThompson et al.l (|2007j ) estimated the missing flux 
component in the J and if— bands from the faint outer 
parts of galaxies and from galaxies below their detection 
limits using a histogram of flux in all pixels associated 
with detected objects. We use a similar method to esti- 
mate the flux missed in the K s — band. In Figure [4] we 
show a histogram of number of pixels versus flux for all 
pixels associated with detected galaxies (object pixels) in 
our MOIRCS K s — band mosaic. Noting the linear trend 
for fluxes ~ 0.005 — 0.4 /xJy, we fit a line to the data over 
this range (blue dashed line). The portion of the his- 
togram used in the fit represents 60% of all object pixels 
in the image and 99% of object pixels with fluxes greater 
than the turnover in the histogr am at ~ 0.005 niy. The 
slope of the linear fit is -0.86. IThompson et al.1 (|2007T ) 
find a slightly steeper value of -1 by simply estimating 
the slope by eye. Assuming that the true flux distribu- 
tion for faint pixels (< 0.005 /xJy) continues along the 
same trend, we extrapolate the linear fit to approximate 
the shape of the histogram when all pixels in the image 
are accounted for. Using this method, we calculate an 
estimate for flux missed in the faint outer parts of galax- 
ies and in galaxies that are below our detection limits. 
We find the missing flux component to be ~ 22% (1.9 
nW m~ 2 sr _1 ) of the total K s — band light from resolved 
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Fig. 4. — Histogram showing Logio of flux in fijy vs. Logio of 
the number of pixels at that flux for all the pixels associated with 
galaxies in the Subaru MOIRCS K s — band image. The blue dashed 
line shows a linear fit from ~ 0.005—0.4 /^Jy, which includes ~ 60% 
of all pixels associated with galaxies and ~ 99% of such pixels 
containing fluxes higher than the peak of the histogram. The slope 
of this line is -0.86. We extrapolate the linear fit toward fainter 
fluxes to the point where all the pixels in the image are accounted 
for (red dashed line). The red hashed area shows the missing flux 
component corresponding to 1.9 nW m~ 2 sr _1 in the K s — band. 



galaxies. 

3.2. The Stellar Component of NIR Light 

The integrated intensity from stars in the NIR (at 
Earth) is comparable to that of resolv ed galaxies (as 
deter mined by models, such as those of iLevenson et all 
120071 or by integrated star counts, as discussed below). 
Thus, in terms of measuring the IGL, it is crucial to be 
able to separate stars from galaxies accurately. In K 10 we 
demonstrate that we are robustly separating stars from 
galaxies in our analysis. Thus, we are able to derive the 
stellar contribution to total NIR light with star counts 
from 15 < JHK S < 24, while groups making unresolved 
measurements of the total NIR EBL must assume stellar 
distribution models to remove the NIR light due to stars 
from their measurements. The integrated starlight from 
our star counts is 14.7 ±2.4 nW m~ 2 sr _1 in J— band, 



band, and 7.6 ±1.8 nW 



10.1 ±1.9 nW m- 2 sr" 1 in H- 
yoT 2 sr _1 in K s — band. 

Four of our five fields are at Galactic latitudes of 
b ~ 50, and one is at b ~ 70 (we average over all 
five in the numbers quoted above). Looking at the 
modeled starlight as a function of Galactic latitude in 
ILevenson et al.l ([2007D , we find the average NIR light due 
to stars in their models at these latitudes is ~ 7 nW m~ 2 
sr _1 in K— band and ~ 16 nW m~ 2 sr _1 in J— band, 
in good agreement with our measured values. Thus, it 
appears that recent total NIR EBL measurements are 
removing the stellar component correctly. However, as 
pointed out in K10, cool dwarf stars may contaminate 
the galaxy counts fainter than JHK ~ 19, suggesting 
that perhaps both IGL measurements and direct total 
NIR EBL measurements retain some level of stellar con- 
tamination, but this would not affect any discrepancy 
between the two. TeV 7— ray observations are concerned 
with the NIR background along the entire path length to 



the blazars and are unaffected by Galactic NIR starlight 
due to its local origin. 

4. COMPARISON WITH PREVIOUS RESULTS 

In Figure [3] we show our IGL results along- 
side a summary of measurements of the NIR IGL 
(black symbols) and EBL (red symbols) over the 
past decade, in cluding results obtained v ia integrated 
galaxy counts (IMadau fc Pozzettl 120001 : iTotani et al.l 
120011 : Thompson 2003) , via total light minus stars and 
zodiacal light (IDwek fc Arendtlll998t IHauser et al.lll998 ; 
Goriian et all 120001: iWrightl 120011: iCambresv et alj |2001 
Matsumoto et al.l 120051 : ILevenson et al.l 12007), and via 



TeV 7— ray opacity (blue hashed area) ()Aharonian et all 
2006). 

The blue hashed area of Figure [5] shows the allowed 
(1 a) NIR EBL intensity (14 ± 4 nW m" 2 sr" 1 ) de- 
rived from High Energy Stere oscopic System ( HESS ) 
observations o f TeV blazars (jAharonian "eTaLl 120061) . 
iMazin fc Raui (|2007T) use 13 TeV blazars and a grid 
of NIR background intensities to further constrain the 
NIR E BL and find approximat e agreement with the re- 
sults of lAharonian et aT . (2006). As noted earlier, how- 
ever, estimates of the NIR EBL from TeV 7— ray opacity 
measurements rely on assumptions about the intrinsic 
SEDs of blazars and that of the EBL, both of which are 
poorly constrained observationally. Our results (black 
diamonds) are some 10 — 20% higher than previous esti- 
mates of the IGL, which puts them closer to EBL esti- 
mates from TeV blazar observations and the most recent 
total NIR EBL measurements. 

A large NIRBE has been found by several groups, 
with perhaps the most striking result being that of 
IMatsumoto et all (|2005f ) . This excess, when combined 
with other NIR and optical background measurements, 
showed an apparent spectral break around 1 /xm. The 
excess was originally attributed to PopIII stars, with 
the break corresponding to the redshifted Lyman limit 
for these stars at z ~ 10. However, a search for the 
possible absorption imprint of this break on the 7— ray 
SEP of blazars di d not find evidence for such a feature 
(IDwek et alJ[2005l). 

IMadau fc SiikT(l2005l) discuss the implications of a 
NIRBE due to massive PopIII stars in the early universe 
and demonstrate that an excess of even a few nW m~ 2 
sr _1 in the J— band would imply energetic requirements 
that are "uncomfortably high" . They note that ~ 5% of 
baryons would be formed into such stars by a redshift of 
9 (roughly twice the fraction converted into stars since 
then) and that the metals created in these PopIII stars 
would have to be hidden in intermediate mass black holes 
(IMBHs) for the universe to avoid exceeding solar metal- 
licity by redshift of 9. They point out that the IMBHs 
themselves coul d produce the NIR ex cess as miniquasars, 
as suggested by ISantos et al.1 (|2002[ ). but that accretion 
onto these IMBHs would dominate the soft X-ray back- 
ground and their mass density in the present day uni- 
verse would exceed that of black holes found in galactic 
nuclei by 3 orders of magnitude. Finally, they show that 
the ionizing flux produced by these PopIII stars would 
exceed by 3 orders of magnitude that required to pro- 
duce the WMAP observed electron scattering depth at 
z = 17. 
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TABLE 2 

NIR Integrated Galaxy Light nW m~ 2 sr _1 ) 



J a H K a h 

This work c 11.7 (+3) ±2.6 11.5 (+3) ±1.5 10.0 (±1.9) ±0.8 

IMadau fc Pozzettil d2000h 9.7^ g 9.0i^7 7.9±i2 

ITotani et an J20M) 10.9 ± 1.1 N/A ' 8.3 ± 0.8 

Thompson (2003) d N/A 7.0 N/A 



a Madau &c Pozzctti (2000) .7— band is centered at 1.1/im while this and other studies 
listed arc centered at 1.25 /^m. 

b This work is done at if a (2.12 /im) while lMadau fc Pozzettil I2000T) and lTotani et al.l 
(2001) are at K(2.2 p.ni). 

c Missing flux estimates are given in parentheses next to our measured values for the 
IGL. The missing flux in the J and H — bands arc taken from Thompson et al. (2007) 
and for K s — band the value is derived as described in Section 13.11 The confidence 
intervals given are the 1 a error estimates from the galaxy counts fits described in 
Scction[3l 
Lower limit. 
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Fig. 5. — NIR IGL and EBL measurements over the past decade as a function of wavelength (f im). This work (black d i amonds) and 
a summary of measure ments of the NIR IGL (other black symbols) via integrated galaxy counts (Madau & Pozzctti 2000; Totani ct al. 
20011; iT hompson 2003 ), total EBL (red symbols) via to t al light minus stars and zodiacal light IIDwek aTA rcndt 1998; Hause r et al .1119981 ; 
Goriian eta l. 2000; Wright 2001; Cambresy ct aTJl2001l ; iMatsumoto et al.|[2(305l ; ILevenson et alll2007l) . and via TeV 7-ray opacity mea- 
surements (blue hashed area) (Aharonian ct al. 2006). The range indicated for the 7— ray work shows the 1 a error range. Note that data 
points at 1.25, 1.6, and 2.2 (im have been shifted slightly in their abscissa values for clarity. Arrows are used to denote upper and lower 
limits and otherwise error bars represent 1 a confidence levels. Our results bring the measurement of the IGL into better agreement with 
TeV 7— ray observations and the most recent total total NIR EBL measurements. The lower limits on our data points show the 1 a error 
estimates associated with the galaxy counts integra tion described in Section [3] while the upper limits show these same 1 a error estimates 
plus the missing flux component derived in Section |3. II and shown in Tabic [2] 
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5. SUMMARY 

We have presented measurements of the IGL in J, H, 
and K s using deep wide-field NLR photometry in com- 
bination with ultradeep MOIRCS K s — band data and 
selections from the literature. These results place the 
best current constraints on the total NIR light from re- 
solved galaxies and serve as a new lower limit to the total 
NIR EBL. While these results are in relative agreement 
with previous measurements, our numbers are 10 — 20% 
higher, bringing them into better agreement with those 
derived from 7— ray experiments and the most recent 
measurements of the total NIR EBL. 

We find the IGL to be roughly isotropic, consistent 
with the expectation of large-scale isotropy in the EBL. 
We confirm that the starlight subtraction for the most 
recent total NIR EBL measurements is correct, so if there 
still exists a foreground subtraction issue in these mea- 
surements, it most likely is associated with the zodiacal 
light. 

While our measurements cannot rule out the existence 
of a NIRBE due to PopIII stars or other exotic early 
universe objects, our new lower limits on the IGL and 
the upper limits foun d from TeV 7— ray experiments 
(|Aharonian et a l. 2006) could now be considered in rough 
agreement with the most recent total NIR EBL mea- 
surements in the J— band, and in near agreement in the 
K— band. 
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